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Abstract 

We use laser-Doppler velocimetry (LDV) experiments and Reynolds-aver aged Navier-Stokes (RANS) simu¬ 
lations to study the characteristic flow patterns downstream of a standardized clockwise swirl disturbance 
generator. After quantifying the impact of the mesh size, we evaluate the potential of various eddy-viscosity 
turbulence models in providing reasonable approximations with respect to the experimental reference. The 
choice of turbulent models reflects current industry practice. Our results suggest that models from the 
k-e family are more accurate in predicting swirling flows than models from the k-uj family. For sufficiently 
resolved meshes, the realizable k-e model provides the most accurate approximation of the velocity magni¬ 
tudes, although it fails to capture small-scale flow structures which are accurately predicted by the standard 
k-e model and the RNG k-e model. Throughout the article, we highlight practical guidance for the choice 
of RANS turbulence models for swirling flow. 
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1. Introduction 

Swirling flows are ubiquitous in many industrial 
applications including furnaces, cyclone separators, 
heat exchangers, and turbines. Yet, efficient nu¬ 
merical prediction with Reynolds-averaged Navier- 
Stokes (RANS) approximations is challenging and 
large-eddy simulations (LES) are still prohibitively 
expensive for practical applications. For strongly 
rotating flows, various standard assumptions used 
in RANS turbulence closures are not expected to 
hold, as discussed, for example, by Jakirlic et 
al. [1]. Consequently, the performance of conven¬ 
tional RANS modeling approaches for predicting 
swirling flow remains elusive and requires further 
verification and validation with experiments. 

Kitoh [2] performed experiments with turbulent 
swirling flow in a straight pipe where the swirl com¬ 
ponent is generated through guide vanes with vari¬ 
able vane angle and a bell-shaped cone at the cen¬ 
ter of the swirl generator. Based on the charac¬ 
teristic tangential velocity distribution, Kitoh [2] 
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categorized the flow field into three regions: wall, 
annular, and core. Due to the streamline curva¬ 
ture, skewed shear directions in the annular region, 
and the resulting anisotropy, Kitoh [2] concludes 
that Reynolds stress models seem more promising 
than eddy-viscosity models in providing accurate 
predictions of swirling flows. This conclusion con¬ 
firms earlier simulation results of Kobayashi and 
Yoda [3] and is also supported by subsequent nu¬ 
merical studies. For example, simulations of Hirai 
et al. [4] show a better performance of Reynolds 
stress models over the standard k-e model and a 
modified k-e model for the prediction of the lami- 
narization phenomenon in swirling pipe-flow. Sim¬ 
ilarly, Cocic et al. [5] report that two-equation tur¬ 
bulence models including RNG k-e, Lauder Sharma 
k-e, and k-uj SST fail to predict experimental data, 
whereas the Launder-Gibson Reynolds stress model 
and the Speziale-Sarkar-Gatski model are found to 
provide better agreement with experiments. 

In contrast, Ghen et al. [6] find that the differ¬ 
ential Reynolds stress model proposed by Launder 
et al. [7] provides quantitatively inaccurate results 
for predicting tangentially injected swirling pipe 
flows, reaching only qualitative agreement with ex- 


Preprint submitted to Elsevier 


October 21, 2015 





perimental results. Further, several workers suc¬ 
cessfully applied turbulence models of the k-e fam¬ 
ily to model swirling flow. For example, Parchen 
and Steenbergen [8] studied turbulent swirling pipe- 
flows for gas-metering applications and reported an 
acceptable behavior of the k-e model, which pro¬ 
vided a better approximation of the experimental 
values compared to an algebraic slip model. Escue 
and Cui [9] found a superior performance of the 
RNG k-e model over a Reynolds stress model for 
moderate swirl intensities in 2D simulations when 
compared to experimental results. Petit et al. m 
successfully applied a k-e model for unsteady sim¬ 
ulations of a swirl generator, achieving good agree¬ 
ment in averaged velocity profiles. In view of the 
available results, a general consensus regarding the 
performance of eddy-viscosity closures for swirling 
flows appears to be lacking. 

1.1. Swirling flows in flow metering 

In this article, we study the swirling flow gener¬ 
ated by a standardized clockwise swirl disturbance 
generator that is used for testing commercial flow 
meters (Figure(a)). Commercial water, heat and 
cooling meters may be exposed to disturbed in¬ 
flows when they are in operation, since the space 
available in realistic installations is often limited 
and prohibits the installation of a straight pipe 
long enough to create fully developed flow condi¬ 
tions upstream of the meter. Common flow dis¬ 
turbances induced by several standard installations 
such as bent pipes or valves may compromise the 
accuracy of meter readings. Consequently, it is nec¬ 
essary to assess the robustness of new products with 
respect to disturbed flow conditions. The stan¬ 
dards EN ISO 4064-2:2014 [U and OIML R 49- 
2:2013 [12] for the type-approval of commercial wa¬ 
ter meters and EN 1434-4:2007 [T3| for the type- 
approval of commercial heat and cooling meters in¬ 
clude standardized tests with artificially generated 
disturbed flows that are meant to emulate disturbed 
flow conditions in realistic installations. One of 
these artificially disturbed flows is generated by a 
standardized clockwise swirl disturbance generator 
(Eigure (a)), which emulates the flow structures 
induced by a double-bent pipe installation. (Eor 
a study of flow structures induced by double-bent 
pipes see, for example, Mattingly and Yeh [13]. ) 

A detailed analysis of flow patterns generated 
by such swirl disturbance generators is not part 
of the standard type-approval procedure, but it is 


a key factor in enabling a meaningful interpreta¬ 
tion of type-approval tests. Eurther, the detailed 
analysis of swirling flows has potential to assist 
the development process of robust meter designs 
and the associated testing and verification facili¬ 
ties. The swirling flow field generated with equip¬ 
ment according to EN ISO 4064-2:2014 [11] and 
OIML R 49-2:2013 [12] is expected to exhibit simi¬ 
lar features as the flow fields in the studies discussed 
in Section Eor example, Eichler and Lederer m 
conducted stereoscopic particle image velocimetry 
(SPIV) measurements downstream of a standard¬ 
ized DN80 swirl disturbance generator, showing 
that the swirling flow is maintained up to 87.OD 
downstream. Similarly, Tawackolian [16] provides 
experimental and numerical results using the k-uj 
model for a DN80 swirl generator and a comparison 
of the flow patterns produced by the swirl generator 
against those after bent pipes. While these studies 
provide valuable insight for large pipe diameters, no 
results appear to be available for smaller pipe diam¬ 
eters. Since the standardized parts for different di¬ 
ameters according to EN ISO 4064-2:2014 [TT] and 
OIML R 49-2:2013 [12] are not exactly self-similar, 
the scalability of available results to smaller diam¬ 
eters remains unclear. 


The present work aims to asses the potential 
of RANS simulations in providing accurate predic¬ 
tions of the flow patterns downstream of a stan¬ 
dardized swirl disturbance generator. In view of the 
discussion in Section available practical guidance 
for choosing turbulence models is ambiguous. To 
elucidate modeling choices for swirling pipe-flow, 
we study different numerical setups and compare 
them against experimental results obtained from 
laser-Doppler velocimetry (LDV) experiments. Our 
choice of turbulence models includes some of the 
most popular mainstream models currently used in 
industry. The aim is to identify best practices to 
achieve realistic and useful simulation results that 
can be used in practical applications such as flow 
meter testing and development. Additionally, we 
discuss our results in the context of available rec¬ 
ommendations regarding the modeling of swirling 
flows. Within the study of different turbulence 
models, we also provide a detailed assessment of 
the impact of the mesh size and the length of the 
computational domain. 
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Figure 1: (a) Standarized swirl disturbance generator ac¬ 
cording to EN ISO 4064-2:2014 [TT] and OIML R 49- 
2:2013 |12| . (b) Measurement grid for the LDV experiments, 
(c) LDV probe and optical access from ILA/Optolution. 


2. Materials and methods 


We use a standardized clockwise swirl dis¬ 
turbance generator according to EN ISO 4064- 
2:2014 [TTj and OIML R 49-2:2013 [12] with an in¬ 
ner diameter oi D = 15.0 mm. The working fluid is 
water and both simulations and experiments are re¬ 
alized with Reynolds number Re = 4.0 • 10^, where 
Re is based on the volumetric velocity w^o\ and the 
pipe diameter such that 


Wyo\D 

Re =-, 

V 


( 1 ) 


and 


_ Q 

'^vol 1 


( 2 ) 


with V the kinematic viscosity of the water, Q the 
volumetric flow rate, and A = the area of the 
pipe cross-section. 


2.1. Experiments 

All experiments are performed on a calibration 
flow bench with an expanded uncertainty of 0.30 % 
(coverage factor k = 2.0) for measurements against 
weight with test-volumes of 3 liter to 100 liter. 
The flow bench is equipped with three magneto- 
inductive (MID) master meters, each one responsi¬ 
ble for controlling the flow in a certain flow-rate in¬ 
terval. The experiments are realized with water at 
temperature T = 20.0 °C, corresponding to a kine¬ 
matic viscosity of z/ = 1004.79 • 10“^ m^/s. During 
measurements, the volume flow-rate Q, the water 


temperature T, and the pressure p are stabilized 
through PID controlled feedback loops. We verify 
the stability of Q, T, and p by logging data from 
the master meters and the corresponding temper¬ 
ature and pressure sensors. Let Qm denote the 
time-averaged master flow-rate and let cfqm de¬ 
note the associated standard deviation providing a 
measure for the stability of the flow-rate. We find 
ctqmIQm ^ d-25 % and that the master meter sig¬ 
nal has the characteristics of random white noise, 
which confirms that there is no preferred timescale 
and no low frequency disturbances that might bias 
the long-time accuracy of LDV measurements. 

We perform measurements with a commercial 
LDV system (ILA/Optolution) in and industrial 
flow laboratory. The flow is seeded with neutrally 
buoyant silver coated hollow glass beads to improve 
the LDV signal. To ensure a fully developed flow 
profile upstream, the swirl generator is located after 
more than lOO.OD of straight pipe and the measure¬ 
ment section is located 12.0D downstream from the 
swirl disturbance generator. For the present study, 
we focus on discussing the flow patterns in the near- 
held until 12.OD downstream from the swirl distur¬ 
bance generator, since this is the how that will en¬ 
ter a how meter in a type-approval test. The mea¬ 
surement grid comprises 241 measurement points 
as shown in Figure [^(b). 

The amount of data acquired at each measure¬ 
ment point is determined by the choice of two ex¬ 
perimental constraints: (I) the maximal number of 
single-point samples rimax and (II) a timeout tmax 
for each single-point measurement on the measure¬ 
ment grid. For the present measurements we choose 

^max 10 and ^max 60.0 s. (3) 

Axial velocity prohles are obtained through a collec¬ 
tion of single-point measurements of the local axial 
mean velocity w over the measurement grid (Fig¬ 
ure [^. To determine nJ, we use the estimator 

1 "" 

(4) 

n 

2=1 

with Wi single-point samples of velocities and n the 
number of samples. The associated dimensionless 
turbulence intensity is 

Tu=^, (5) 

w 

where 

/ n \ 1/2 

(Jw = ( - wf j (6) 
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is the standard deviation of samples wi. To esti¬ 
mate the reliability of the estimator at each spa¬ 
tial measurement point, we determine the associ¬ 
ated standard error 


is the turbulent kinetic energy. In view of pTj 
and the RANS equations with the gradient- 
diffusion hypothesis and the eddy-viscosity hypoth¬ 
esis are 


(yw = (ywl^fn. {7) 

The standard error provides an uncertainty es¬ 
timate for the estimator of the mean velocity 
The mean velocity components u in x direction 
and V in y direction are determined analogously. 
To obtain all velocity components, three consecu¬ 
tive LDV measurements are required. Notice that 
the velocity vector u = (H, U, wY is the experimen¬ 
tal equivalent of the Reynolds-averaged velocity de¬ 
fined in Section [221 


2.2. Simulations 

The numerical simulations are performed using 
a Reynolds-averaged Navier-Stokes (RANS) ap¬ 
proach. Let a superposed bar denote Reynolds- 
averaging, so that (see, for example Pope m) 

u = u + u', (8) 

where u is the RANS velocity and u' is the fluctuat¬ 
ing velocity satisfying u' = 0 . Applying Reynolds- 
averaging to the incompressible Navier-Stokes (NS) 
equations yields the RANS equations 

u = — ^gradp + z/Au — divr, 
divu = 0, 

where a superposed dot denotes material time- 
differentiation following u, p denotes the averaged 
pressure, p denotes the fluid density, and 

T = u' (g) u' (10) 



is the Reynolds stress tensor. The closure prob¬ 
lem in RANS modeling amounts to finding expres¬ 
sions for without using the fluctuating velocity 
u'. The gradient-diffusion and the eddy-viscosity 


hypotheses assume that the deviatoric part of 10 


is proportional to the averaged stretching tensor 


I (grad u + (gradu)^) and yield the closure 


D 


approximation 


- T + = J/tD, (11) 

where z/t is the eddy-viscosity, I is the identity ten¬ 
sor, and 


k = 


u 


/|2 


( 12 ) 


u = (z/ + z/t)Au - ^grad(p -h ^pk), 
divu = 0. 

Notice that k can be absorbed into a modified av¬ 
eraged pressure. 

In this article, we consider various two-equation 
eddy-viscosity models (Table that determine z^t 
through solving two additional model transport 
equations. The first additional model transport 
equation is for the turbulent kinetic energy k and 
the second additional model transport equation de¬ 
pends on the model family. For models from the k-e 
family, the second additional transport equation is 
for the turbulent dissipation 

e = 2z/D' : D' (14) 

where D = | (grad u'+ (grad u')^) is the fluctuation 
stretching tensor. Similarly, for models from the k- 
uj family, the second additional transport equation 
is for the turbulent frequency uj = e/k. Despite 
the vast selection of available turbulence models, we 
focus on popular mainstream models that are com¬ 
monly used in industrial application. More exotic 
models specifically developed for swirling applica¬ 
tions might well yield improvements in performance 
with the downside of giving up universality, which 
can be problematic in industrial applications. For 
a discussion of model performance, also see Section 
5. 

2.2.1. Numerical solution 

We use the open-source CFD code Open- 
FOAM [18] and assess five different two-equation 
eddy-viscosity models, as summarized in Table 
Both convective and diffusive terms are discretized 
with a second order Gaussian linear scheme. The 
near-wall region is modeled with enhanced wall 
functions available in the standard OpenFOAM dis¬ 
tribution. A SIMPLE algorithm is used for the cou¬ 
pling between the pressure and the velocity equa¬ 
tions. The pressure equation is solved with a geo¬ 
metric agglomerated algebraic multigrid (GAMG) 
method, and all other equations are solved with the 
smoothSolver solver of OpenFOAM using a Gauss- 
Seidel smoother. 

The computational domain has a total length 
of 400.0 mm with the swirl disturbance generator 
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Table 1: Turbulence models and references of the specific 
model formulations used in the implementation in Open- 
FOAM._ 


Model 

Main reference(s) 

standard k-e 

Launder and Spalding |19| 

RNG k-e 

Yakhot et al. |20| 

realizable k-e 

Shih et al. |21| 

k-uj 

Wilcox [22] 

k-u SST 

Menter and Esch |23| & Hellsten |24| 


placed at the center. Consequently, the distance 
from both the inlet and outlet to the swirl generator 
is approximately 13.0D. We use hex-dominant un¬ 
structured meshes generated with the OpenFOAM 
meshing tool snappyHexMesh. All meshes are gen¬ 
erated with the mesh quality constraints summa¬ 
rized in Table |B3] in [Appendix B| 


(a) (b) 




-©-standard k-e 
-b-RNG k-e 
—I— realizable k-e 
->i-k-uj 
-^k-io SST 
—Theory Moody 


Figure 4: (a) Pressure along the centerline for the different 
turbulence models, (b) Darcy friction factor for the pipe 
section between —3.0D and —l.OD. The inlet is located 
at x/D ^ —13.0, the swirl disturbance generator is located 
at x/D = 0.0 (shaded area), and the outlet is located at 
x/D ^ 13.0. 


3. Results 

3.1. Impact of the mesh size and domain length 

To identify an optimal trade-off between simu¬ 
lation time and accuracy of results, we assess the 
impact of the mesh size and the impact of the length 
of the computational domain on the solution of dif¬ 
ferent turbulence models. To that end, we run all 
turbulence models with five different mesh sizes. 
To quantify the effect of the domain length, we 
first perform simulations with a 20.OD inlet and 
outlet corresponding to a total domain length of 
800.0 mm. In Figure we compare results of both 
inlet lengths for three selected turbulence models. 
Panel (a) of Figure shows the axial component w 
of the velocity 12.OD downstream from the swirl 
disturbance generator. Further we visualize the 
magnitude of the secondary flow 

Vxy = \/m^ +V^, (15) 

as shown in panel (b) of Figure Visual com¬ 
parison of both results shows that the character¬ 
istic flow patterns are maintained across different 
domain lengths although there are small-scale dif¬ 
ferences in the velocity field. To assess these differ¬ 
ences in more detail, we compare velocity profiles 
over r/R ol each case upstream and downstream 
from the swirl generator. We compute averages of 
10 individual linear profiles with equidistant angu¬ 
lar spacing of 18° in analogy with the experimental 
measurement grid shown in Figure [^(b). The com¬ 
parison of the averaged axial profiles upstream from 
the swirl generator shows that both RNG k-e and 


k-uj SST models provide flow profiles closer to the 
fully developed reference with the longer inlet, while 
the realizable k-e model appears to provide similar 
results for both inlet lengths (Figure]^ (c)). In the 
cross-section 12.OD downstream from the distur¬ 
bance generator, the RNG k-e model shows a slight 
attenuation of the axial peak velocity for the long 
inlet and the k-uj SST model exhibits a smoother 
profile for the long inlet. The corresponding Vxy 
profiles are shown in panel (e) of Figure Both 
the realizable k-e model and the RNG k-e model 
show almost identical Vxy for the long and the short 
computational domains. For different lengths of 
the computational domain, the k-uj SST model pre¬ 
serves main features including peaks and gradients 
close to the wall and the pipe center, but shows 
small differences in the Vxy profile for intermedi¬ 
ate r/R. These results suggest that the impact of 
the domain length in the considered cases is minor. 
Further, this comparison illustrates, that the profile 
upstream from the swirl disturbance generator has 
a minor influence on the flow patterns downstream, 
suggesting that upstream disturbances do not in¬ 
fluence the characteristic patterns of the swirling 
flow field downstream. Therefore, the study is per¬ 
formed using the shorter computational domain, 
which allows a reduction of the number of cells and 
the computational time. 

Next, we study the influence of the mesh-size on 
the solution. Figurej^shows the secondary flow pre¬ 
dicted by different turbulence models and different 
meshes. The standard k-e model and the RNG k- 
e model show little dependency on the mesh size, 
maintaining similar flow patterns across different 
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(e) 


1 -1 


0 

vjR 


— RNG k-e short 
---RNG k-6 long 
— realizable k-e short 
— realizable k-e long 
k-LU SST short 
k-uj SST long 



Figure 2: The length of the computational domain has minor influence on the downstream solution: Comparison of uJ/iCvoi and 
^xy/w^oi for the long and short computational domains for three different turbulence models, (a) Contour plots oiw/w^oi at 
12.OD downstream, (b) Contour plots of Vxy/w^ol at 12.OD downstream, (c) Averaged w/w^oi profile at 2.011 upstream, (d) 
Averaged w/w^oi profile at 12.011 downstream, (e) Averaged Vxy/wvo\ profile at 12.011 downstream. Note that the secondary 
flow upstream from the swirl generator is zero. 


mesh sizes. In contrast, the realizable k-e model ex¬ 
hibits a significant difference between solutions on 
coarser and finer meshes. While the coarser meshes 
show small-scale flow patterns, the finer meshes pre¬ 
dict a homogeneous swirling flow without charac¬ 
teristic small-scale patterns. In general, the models 
from the k-e family tend to predict smaller mag¬ 
nitudes of Vxy with increasing mesh sizes, whereas 
models from the k-uj family tend to predict higher 
magnitudes of Vxy with increasing mesh size. Un¬ 
der the assumption of monotonic convergence, this 
suggests that swirling flows computed with mod¬ 
els from the k-e family should be interpreted as a 
lower estimate for the secondary flow, whereas re¬ 
sults computed with models from the k-uj family 
should be interpreted as an upper estimate for the 
secondary flow. 

Further, we study the convergence of the pressure 
drop with increasing mesh-size. The pressure along 
the center-line of the pipe is shown in Figure]^ (a) 
for all considered turbulence models. As expected. 


there is a significant pressure drop across the swirl 
disturbance generator, while the pressure loss is ap¬ 
proximately linear with respect to the ^-coordinate 
for the pipe sections before and after the swirl dis¬ 
turbance generator. To compare the pressure loss 
in the straight pipe with theoretical references, we 
compute the Darcy-Weisbach friction factor 


/ = 


2D 


Ap, 


(16) 


where p is the density of water, L is the associated 
length of the pipe-section (from —3.01) to — l.OD, 
to consider the section where the flow is most de¬ 
veloped), and Ap is the associated pressure loss. 
The friction factor is computed for all turbulence 
models and compared to the theoretical reference 
obtained from the Moody diagram [25]. Figure]^ 
(b) shows that the pressure drop over the straight 
pipe downstream from the swirl disturbance gener¬ 
ator is overpredicted by all considered turbulence 
models when compared to the theoretical reference 
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Figure 3: Mesh dependency of different turbulence models: Comparison of Vxy/w^^o\ at 12.0D predicted by various models and 
mesh sizes. 


from the Moody Diagram. The standrad k-e model 
and the RNG k-e model appear to converge to a 
value that is higher than the theoretical reference 
as the mesh size increases. Moreover the realizable 
k-e model shows a significant overprediction, where 
values for the finer meshes do not agree with the 
theory nor with the other models. Conversely, mod¬ 
els from the k-uj family show a very close agreement 
with the theory and appear to converge towards the 
theoretical value for increasing mesh size. For the 
discussion in the remainder of the article, we use 
the results obtained with a mesh size of 20 • 10^ 
cells. 


3.2. Velocity fields 

For the comparison of flow patterns at different 
cross-sections downstream of the swirl disturbance 
generator, we focus on three of the models since the 
structures predicted by the RNG k-e^ the standard 
k-e models, and both k-oo models are found to be 
qualitatively similar (see, for example. Figure [^. 
The flow patterns of the secondary flow predicted 
at four locations are shown in Figure and iso¬ 
surfaces corresponding to w = Wyoi of the domain 
downstream of the swirl disturbance generator are 
shown in Figure 

All considered models predict eight characteristic 
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Figure 6: Isosurfaces oiw = Wyoi representing the evolution 
of the flow downstream of the swirl generator for all the 
considered turbulence models. 


flow patterns at the cross-section 2.0D downstream, 
which are associated with the eight blades of the 
swirl disturbance generator. At 6.01), the realizable 
k-e model does not show any small-scale flow pat¬ 
terns in Vxy^ but a radially homogeneous swirling 
flow. In contrast, the k-uj SST model and the k-e 
model still predict eight patterns with elevated 
The homogeneous swirling flow of the realizable k- 
e model is maintained further downstream with a 
continuous reduction of Vxy. For the k-e model, 
the initial structures start dissipating and merging 
at around lO.OD, resulting in only four visible flow 
structures at the 12.0D measurement cross-section. 
The k-uj SST model maintains the eight character¬ 
istic flow patterns throughout the entire simulation 
domain, although at the 12.0D measurement sec¬ 
tion the dissipation of four of them starts to become 
visible. 

Next, we compare the simulations with the ex¬ 
perimental results at the section 12.0D downstream 
of the swirl disturbance generator. Figure shows 
the secondary flow obtained from simulations along 
with the corresponding experimental results. Vi¬ 
sual inspection of the contour plots suggests that 
the standard k-e model and the RNG k-e model pro¬ 
vide the best approximations of the experimentally 
observed four characteristic patterns in the outer 
region of the pipe. Conversely, the realizable k-e 
model provides a more realistic approximation of 
the velocity magnitude close to the center of the 
pipe. 

For the comparison of individual profile paths, we 
use 10 averaged profiles as described in Section [tT| 
The individual measured velocity profiles are shown 


in Figure along with the corresponding standard 
error. The standard errors associated with the aver¬ 
aged profiles are computed through error propaga¬ 
tion calculus using the standard error of individ¬ 
ual profiles. For the averaged axial velocity profile, 
we compute the associated standard error 


/ m \ 


(17) 


where m = 10 is the number of individual profiles. 
Similarly, for the secondary flow, assuming uncor¬ 
related au and we find 


^ xy — 


1 

m 



u‘f + v‘f 


1/2 


(18) 


The averaged profiles are shown in Figure The 
shaded area indicates the associated standard er¬ 
ror of the experimental results with a coverage fac¬ 
tor K = 2.0. The profiles of w indicate that the 
realizable k-e model is the only model that does 
not predict a spurious pronounced velocity peak at 
the pipe center, but an almost flat plateau-like pro¬ 
file that closely approximates the experimental data 
(Figure ID (a)). 

Similarly, the profiles of Vxy show that the near- 
held how produced by the standardized swirl dis¬ 
turbance generator is similar to solid-body rotation 
(Figure (b)). Importantly, the how held gen¬ 
erated by the standardized swirl disturbance gen¬ 
erator according to EN ISO 4064-2:2014 m and 
OIML R 49-2:2013 [12] appears to lack a transition 
of the how from solid body rotation to Rankine- 
vortex type how as reported by Vaidya et al. [26] 
for a swirling how generated by a rotating honey¬ 
comb. All turbulence models appear to overpredict 
Vxy and, as previously seen in Figure the high¬ 
est values are given by the models of the k-uj fam¬ 
ily. Overall, the best prediction of the experimental 
data is provided by the realizable k-e model which 
appears to capture the mean how velocity more ac¬ 
curately than the other considered models. 


3.3. Performance indicators 

Performance indicators are useful integral met¬ 
rics to quantify how conditions. Following Yeh and 
Mattingly m and Muller and Dues |28], we com¬ 
pute the swirl angle 0, the prohle factor i^p, the 
asymmetry factor and the turbulence factor 
Ktu- a detailed dehnition of the performance indi¬ 
cators is provided in [Appendix A[ 


8 














Figure 5: Evolution of the secondary flow downstream of the swirl generator: Contour plots of Vxy/w^oi for various cross-sections 
downstream. 


(a) experimental 




standard k-e 



Figure 7: Comparison of simulations and experiments 12.OD downstream from the disturbance swirl generator: Vxy/w^oi from 
LDV experiments (a) and predicted by the different turbulence models (b)-(f). 





In Table we compare performance indicators 
from computational and experimental results as 
well as available experimental results of Eichler and 
Lederer [15] with a DN80 pipe diameter and a flow 
rate of 96.221 m^/h (Re 4.0 • 10^). The values of 
performance indicators are reported as mean values 
of 10 profiles shown in Figure (b) along with the 
associated standard deviation. According to the 
definition A.4, the asymmetry factor is signed. 


However, to report a meaningful non-zero mean, we 
average the magnitude of Kg,. For the simulations, 
the turbulence factor K^u is computed using the 
turbulent kinetic energy k that is part of the so¬ 
lution of the corresponding turbulence model, such 
that 


Tu = 


VW3 


(19) 
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Figure 8: Experimental velocity profiles at 12.OD downstream for the three Cartesian components with standard errors (shaded 
area): (a) ^-component w/w^o\, (b) y-component F/iCvob and (c) a:-component n/iCvol- 
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— standard k-e 
— RNG k-s 
— realizable k-s 
— k-uj 
— k-uj SST 


Figure 9: Comparison of simulations and experiments 12.0 i9 downstream from the swirl disturbance generator: (a) Ave rag ed 
profiles ofw/w^oi profile and associated standard error |17| (b) Averaged Vxy/wyoi profile and associated standard error |18| 


Notice that corresponds to a three-dimensional 
turbulence factor, while the Ktu value of the ex¬ 
periments corresponds to individual axial profiles. 


The values in Table [2] indicate that the realizable 
k-e model provides low standard deviations in per¬ 
formance indicators. This is consistent with the cir¬ 
cular symmetry of the predicted flow fields, which 
results in almost identical individual profiles, as dis¬ 
cussed in Section 3.2 The profile factor of the 
realizable k-e model shows good agreement with the 
LDV results, while all other models predict a higher 
value for this parameter. This substantiates the re¬ 
sults discussed in Section |3.2[ where the realizable 
k-e model is found to provide the best prediction of 
the experimental axial profile. 


The asymmetry factor iCa obtained from the sim¬ 
ulations is significantly lower than the experimen¬ 
tal value. Consequently, none of the models shows 
potential to provide an accurate estimation of this 
parameter. However, the experimental setup is al¬ 
ways subject to small imperfections, which results 
in non-symmetric flow profiles, whereas the com¬ 


putational model is fully symmetric by definition. 
Further, previous investigations (see, for example, 
Eichler and Lederer m) reported that the swirling 
flow field may develop instabilities resulting in sym¬ 
metry breaking at a certain distance downstream. 
The present computational model does not capture 
the effect of symmetry breaking. The stability of 
the swirling flow could be assessed computation¬ 
ally through small perturbations, for example in the 
boundary conditions. However, such an investiga¬ 
tion goes beyond the scope of this article and should 
be investigated in future studies. 

The turbulence factor K^u is approximated rea¬ 
sonably well by the considered turbulence models. 
The best prediction is provided by the RNG k- 
e model and the standard k-e model. The val¬ 
ues obatined from the k-oo model and the k-uj SST 
model are lower than the experimental values. The 
realizable k-e model appears to fail in capturing the 
behaviour of giving a considerable overpredic¬ 
tion. 

The swirl angle is well predicted by the three k- 
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e models, whereas the k-u models predict values 
that exceed those obtained from the LDV experi¬ 
ments. In the present scenario, the realizable k-e 
model gives the closest value to the experiments, 
which supports the accuracy of this model in pre¬ 
dicting the experimental radial profiles as discussed 
in Section |3.2[ In general, simulations with turbu¬ 
lence models from the k-e family appear to provide 
the most reasonable predictions of the performance 
indicators. Yet, none of the considered models gives 
an accurate estimation of all four parameters. 

Comparing our experimental and computational 
results with the experimental results of Eichler and 
Lederer [15], we see a reasonably close agreement 
for iCp, i^Tu: 0- The observed discrepancy in 

i^a can be influenced by the experimental setup of 
the LDV system where the pipe center needs to be 
determined by the operator. This procedure is not 
very robust and can result in small displacements 
from the actual center causing asymmetric flow pro¬ 
files. Further, the conditions for the swirling flow 
field to become unstable remain elusive, as dis¬ 
cussed by Tawackolian [16] and Bilker [29]. With¬ 
out additional information, the distinction between 
actual asymmetries in the flow field and asymme¬ 
tries due to the uncertainties associated with the 
determination of the pipe center is inaccessible. 

4. Discussion of model performance and 
modeling assumptions 

In conclusion, eddy-viscosity models only provide 
limited accuracy for reliable predictions of swirling 
flow. However, these models are still dominant in 
industrial applications since other methods such as 
large eddy simulations (LES) still have a significant 
computational cost, which is presently not achiev¬ 
able in the industrial time-scale. 

To discuss how different modeling assumptions 
of the considered eddy-viscosity models influence 
the predictions of swirling pipe-flow, we take the 
standard k-e model as a baseline reference. 

The RNG k-e model includes an extra term in the 
e equation that is expected to enhance its perfor¬ 
mance in flows with high strain rate and streamline 
curvature (see, for example, Escue and Cui [9]). Al¬ 
though there is a considerable streamline curvature 
in the present scenario, we do not detect major im¬ 
provements with respect to the standard k-e model. 
However, a potential improvement of the RNG k-e 
model by using an empirical modification of the tur¬ 
bulent viscosity is part of various commercial codes 


and was applied by Escue and Gui [9] and Saqr et 
al. [30]. The implementation of such a modification 
in OpenFOAM might well yield further improve¬ 
ments, but goes beyond the scope of this article. 

Similarly, the realizable k-e model aims to pro¬ 
vide improvements with respect to the standard k-e 
model in swirling flows through taking into consid¬ 
eration the effects of mean rotation on computing 
the turbulent viscosity. In view of this property, 
the ability of the realizable k-e model to capture 
the velocity magnitudes more accurately than the 
standard k-e model appears to be consistent with 
the modeling assumptions. 

The k-uj model by Wilcox [22] is designed to im¬ 
prove the performance of the standard k-e model 
in the viscous layer of the near-wall regions and 
in adverse pressure gradients, promising more ac¬ 
curate results for free shear flows and separated 
flows. Additionally, the k-uj SST model combines 
the k-uj model and the standard k-e model through 
a blending function that prompts the use of the k-uj 
model in the near wall regions and the k-e model 
in regions far from the wall. The k-uj SST model 
was originally formulated by Menter m with the 
goal of modeling aerodynamics flows. Menter and 
Esch [23] presented a modification of the k-oo SST 
model which includes a new definition of the eddy 
viscosity using the strain rate instead of the vor- 
ticity. This modification is aimed at extending the 
applicability of the model to a wider range of flows. 
Yet, in the present scenario, we see no improve¬ 
ments on using the k-uj model or the k-uj SST model 
rather than the standard k-e model. With respect 
to the experimental results, both models give the 
most unsatisfactory results for predicting realistic 
velocity profiles in the 12.OD measurement plane. 
Gonversely, we find that the k-uj model and the k- 
uj SST model provide realistic values of the pressure 
loss for sufficiently fine meshes. Further, the mod¬ 
eling assumptions suggest that the k-uj models are 
likely to provide a realistic prediction of the flow in 
the pipe section close the swirl disturbance genera¬ 
tor and inside the swirl disturbance generator, were 
the influence of the wall region becomes more rele¬ 
vant. Yet, these modeling assumptions render the 
k-uj models deficient for predicting the evolution of 
the flow further downstream, where the influence of 
the wall is less dominant. 
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Table 2: Performance indicators at the cross-section 12.OD downstream from the swirl disturbance generator. 



KA-\ 

iCa[%] 

^Tu[-] 

b[deg] 

standard k-e 

1.173 ±0.131 

0.062 ± 0.035 

1.792 ±0.004 

12.269 ±0.550 

RNG k-e 

1.186 ±0.174 

0.051 ± 0.045 

1.910 ±0.007 

12.514 ±0.856 

realizable k-e 

0.438 ± 0.005 

0.002 ±0.001 

3.867 ±0.000 

11.292 ±0.017 

k-(jj 

1.163 ±0.152 

0.241 ±0.325 

1.176 ±0.004 

15.068 ±1.062 

k-u SST 

1.360 ±0.207 

0.309 ±0.356 

1.340 ±0.006 

15.293 ±1.521 

Present experiments 

0.470 ±0.196 

0.975 ±1.120 

2.071 ±0.054 

11.426 ±0.561 

Experiments 13.OD down¬ 
stream with DN80 and 
S3 4 ■ 10^ [H] 

0.61 

0.26 

2.14 

13.0 


5. Conclusions 

We assessed the potential of five different turbu¬ 
lence models to predict flow patterns downstream of 
a standardized swirl disturbance generator with dif¬ 
ferent meshes and mesh topologies. The numerical 
results were validated through a systematic com¬ 
parison with LDV experiments. For sufficiently fine 
meshes, the realizable k-e model gives the most rea¬ 
sonable prediction of both axial and in-plane veloc¬ 
ity profiles, when compared to experimental results. 
However, the realizable k-e model fails to predict 
the characteristic small-scale flow patterns deter¬ 
mined experimentally, which are well-captured by 
the standard k-e model and the RNG k-e model. 
The pressure drop over the straight pipe-section be¬ 
fore the swirl generator is overpredicted by models 
from the k-e family when compared to a theoreti¬ 
cal reference determined from the Moody diagram. 
In particular the realizable k-e model appears to 
provide elevated values for the friction factors. On 
the other hand, the k-uj models predict values that 
are closest to the theoretical reference. We con¬ 
clude that, in general, k-e models appear to pro¬ 
vide more accurate predictions of the velocity fields 
but fail to provide realistic predictions of the pres¬ 
sure drop. Conversely, k-uj models give a better 
prediction for the pressure drop. Overall, the per¬ 
formance of the investigated eddy-viscosity models 
appears to be limited, which is in agreement with 
expectations (see, for example, Jakirlic et al. [T]). 
However, the usage of these models still reflects cur¬ 
rent industry practice and this detailed study in¬ 
cluding experimental validation facilitates a more 
targeted model selection and guidance for the in¬ 
terpretation of results. Additionally, a comparison 
of the present results with other references with dif¬ 
ferent pipe diameters and flow rates shows that the 
performance indicators exhibit similar values. This 
suggests that, despite small differences in the swirl 


disturber geometries, the How fields generated at 
different pipe diameters and How rates are approxi¬ 
mately self-similar. Consequently, small differences 
in the geometry of the swirl disturbance generators 
have only secondary effects on the How field. 


Appendix A. Definition of performance in¬ 
dicators 


Appendix A.l. Profile factor 

Following Yeh and Mattingly [27] , the dimension¬ 
less profile factor Kp is defined as 


d^p,m 


with 


(A.l) 


R 


d^p,m — 

w^o\D J 

-R 


and 


K, 


p,s 


It 

J(^rn,s-Ws)dr, 


(A.2) 


(A.3) 


where = w{r = 0) is the velocity at the pipe 
center, Wm,s is the velocity of the norm profile at 
the pipe center, and Wg is the velocity of the norm 
profile. The profile factor is a measure for peakness 
{Kp > 1) or flatness {Kp < 1) of measured velocity 
profiles with respect to standard profiles such as 
Hagen-Pouiseuille for laminar flow or Cersten and 
Herwig jSH |33] for turbulent flow. 


Appendix A.2. Asymmetry factor 

Following Yeh and Mattingly [27] , the asymmetry 
factor 


K, 


R 

J rw dr 
1 -R 

D R _ 

J w dr 

-R 


(A.4) 
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quantifies the relative radial displacement of the 
center of gravity of the area under the flow profile 
with respect to the pipe center. 


Appendix A.3. Turbulence factor 

Each LDV point measurement is a collection of 
a large number of bursts resulting in a histogram 
(or probability density function) for the axial ve¬ 
locity component. The level of dispersion (i.e. 
the standard deviation) of this histogram quanti¬ 
fies the turbulence intensity As discussed by 
Durst et al. [34] and generalized by Pashtrapan- 
ska [35] , the turbulence intensity in the core re¬ 
gion —0.2 < r/R < 0.2 can be estimated as 

/ \ 

Tu, = 0.13 ( ) (A.5) 

V W^vol J 

for > 4500. The turbulence factor Ktu is 

defined as 


i^Tu 


TUniax 

Tu, 


(A.6) 


where Tumax is the maximum of[^in the core region 
-0.2 < r/R <0.2. 


Appendix A.4- Swirl angle 

Following Yeh and Mattingly [36], the level of 
swirl can be measured quantitatively through the 
maximal swirl angle 

(/. = arctanf<A^). (A.7) 

V '^vol J 

However, the precise debt ion of the swirl angle may 
vary slightly depending on the author. Geometri¬ 
cally, the swirl angle |A.7| is the angle between the 
ideal velocity vector and the actual velocity vector 
with swirl. 
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